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r^, ' Abstract 

^vq ! We present analytic methods for extracting a class of bulk geometries given infor- 
^ ' mation of certain physical quantities in the boundary CFT. More specifically we look 
J£j^ ■ at singular correlators and entanglement entropy in the CFT to provide information 
\Q ', of null and spacelike geodesies repectively in the bulk. We show that static spher- 
ed I ically symmetric, asymptotically AdS spacetimes which do not admit null circular 
t~^ ■ orbits can be fully recovered, and that any spacetime can be recovered up to the local 
^^ , maximum of the potential. We provide analytical and numerical examples to verify 
Q ', the methods used. 
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1 Introduction 

The AdS/CFT correspondence relates two objects: 1) Type IIB superstring theory on 
AdSs X S^ and 2) A/" = 4 SU(N) super- Yang- Mills theory m 3 + 1 dimensions [1, 2, 3]. 
The AdS/CFT correspondence is holographic, since it relates theories living in different 
dimensions. One consequence is that a d + 1 dimensional asymptotically AdS gravitational 
theory is "dual" to a d dimensional field theory which lives on the boundary of this space- 
time. For example, it is well known that a spacetime containing a large black hole has a 
thermal CFT dual, where its temperature is the Hawking temperature [4]. 

To gain a deeper understanding of many issues that arise in quantum gravity, it has been 
useful to probe the AdS/CFT correspondence. Many problems of qauntum gravity manifest 
themselves in what happens beyond the horizon of black holes. Much research has been 
done into extracting information from inside the horizon of black holes from knowledge of 
the boundary field theory correlators [5, 6, 7, 8, 9, 10, 11]. This has motivated original work 
by V. Hubeny et.al.[12], to investigate singularities in the boundary correlation functions. 
They discovered that a boundary correlation function G{x,x') = {0{x)0{x')) is singular 
if there exists null geodesies connecting x and x'. The singularities arising from such a 
relation are known as "bulk-cone singularities". This relationship was then used to answer 
questions such as: can we see horizon formation directly in the gauge theoryl 

Complementary to investigating local properties of the CFT, one can also study non- 
local quantities. Then one can ask questions about what physical quantities would they 
correspond to in the bulk theory . One such quantity is the entanglement entropy Sa- In 
particular, it has beed recently shown [13, 14] that the entropy Sa in a d+ 1 dimensional 
CFT can he determined from the d dimensional minimal surface 7a (in the bulk AdS) whose 
boundary is given by the d — 1 dimensional manifold d'~fA = dA. The entropy is given by 
applying the usual area/entropy relation. In the case of AdSs, the minimal surface is 
one dimensional, and is in fact a spacelike geodesic with endpoints on the boundary. The 
entropy in this case is just derived from the "area" of the spacelike geodesic, i.e. the proper 
length C 

Using the relation shown in [12], numerical work was done to reconstruct certain static, 
spherically symmetric bulk metrics given the location of bulk-cone singularities [15, 16]. 
This paper complements the previous work by attempting to reconstruct the bulk metrics 
using inversion techniques and thus provide an exact solution. This is done by noting that 
the location of the boundary singularities can be expressed in terms of two functions, A(j){a) 
and At{a), where if G{x, x') is some singular correlator between the points x = (t, Q) and 
x' = {t',Q'), then t — t' = At and Q — Q' cc A0 (see figure. 1). These two functions can 
in turn be completely determined by the bulk metric, thus an inversion is required to 
determine the bulk metric. 



In this paper we will consider two types of spacetime. In section 2 we will descibe static 
spherically symmetric spacetimes described by one function /(r). We will then introduce 
the inversion technique following from the ideas set out in the introduction. Analytical 
and numerical examples are then given to test the validity of the inversion. In section 3 we 
generalise to static spherically symmetric spacetimes decribed by two functions /(r) and 
h{r). It will be noted that information of the endpoints of null geodesies will not be enough 
to fully determine the metric. However, supplementing the bulk-cone singularity data with 
the entanglement entropy for arbitrary CFT regions will turn out to suffice to extract the 
bulk metric. Since entanglement entropy in the two dimensional CFT can determine the 
proper length C of spacelike geodesies, and the proper length is expressed in terms of the 
metric, it will be shown that this is enough information in the specific case of an AdSa 
bulk to completely determine the metric exactly. Analytical and numerical examples are 
again given to test the validity of, in this case, two inversions. In section 4 we will give 
a discussion of the results. It will be noted that although we have extracted spherically 
symmetric asymptotically AdS spacetimes only in the case of 3 dimensions, it is completely 
valid for any dimensions. We will also discuss limitations of the inversion, and future areas 
that can be explored. 
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Figure 1: (a) is a 2 + l-dimensional plot of boundary to boundary null geodesies in an AdS- 
black hole geometry with metric ds"^ = —f{r) dt^ + jr^r + r'^dVt\ and /(r) = 1 + r^ — j^-, 
where n — 1 angular directions have been projected out. To plot the null trajectories, 
the spacetime has been compactified using the coordinate transformation r = tan(f), thus 
the boundary is a timelike cylinder with radius f = |. (6) illustrates the quantity At (a) 
via projection of the null geodesies on the r-t plane, where a = -j parameterises the 
null curves, (c) illustrates the quantity A0(a) via projection of the null geodesies on the 
r-(/) plane. Notice that the boundary null curve (a = 1), we have At(Q:) = A(j){a) = 
TT, which is as expected from a pure AdS calculation. (Colour coding for values of a: 
{1, 1.3, 1.7, 1.85, 1.87} = {Black, Red, Green, Blue, Purple}) 



2 Static spherically symmetric spacetimes decribed 
by one funcion f{r] 



An important class of static spherically symmetric spacetimes in r?, + 2 dimensions have 
metrics, in global coordinates, of the form 

ds' = -f{r)dt' + ^^+r'dnl (1) 

Putting the condition that as r ^ oo, /(r) ^ r^ + 1 restricts the spacetime to be asymp- 
totically AdS. If we set /(r) = 1 + r^ — p{r), then for the spacetime to be physically 
meaningful, p{r) must be 0(r") where n ^ —2. 

Fixing one endpoint for all null geodesies on the boundary, calculation of the angular, 
A(f){a), and time, At{a) (see figure. 1), separation of null geodesies (labelled by the ratio 
a = j), with endpoints on the boundary (see Appendix A) yeilds the two expressions 

Atia) = 2 / , (2) 

and 

1"^ dr 

A0(a) = 2 / , (3) 



r„^^„{a) ^2 /q,2 _ /W 



where rmin{ct) is the minimum radius reached in the bulk by a null geodesic, and is defined 

Examining the two terms, there is a relation between A(f){a) and At{a) (see Appendix 
B for derivation) which goes as, 

a^At{a) = ^A(p{a) (4) 

da da 

Although it appears that we have two functions worth of data to extract f{r), they are 
not independent from the equation above. Thus we will use A0(a) only for the extraction. 

2.1 Extracting /(r) 

Examining the definition for A(p{a), it is first noted that this is an integral equation of the 
first kind. The simplest equation of this form is a{x) = J^ b{r) dr, which has the trivial 
solution b{r) = a'{r). The extraction of /(r) is complicated by the fact that equation (3) is 
non-linear in /(r). Non-linearity in any type of equation makes finding an analytic solution 
much harder, so it would be useful to find a quantity which is linear in the integrand and 
invert that instead. This can be done by identifying linear equations which control the 



trajectory of null geodesies. One such equation is the "energy equation". From Appendix 
A we have the equation for f, 



.2 2 fir 



P = a^-^ (5) 

/(r) 
By inspection, we can identify the potential V{r) = — —, to get, 

f2 + v{r) = a^ (6) 

Instead of integrating over the quantity r, let us perform a substitution to integrate 
over the quantity V{r) instead. If we let rmj„(a) = x, we see an immediate simplification 
for the expression of A(p{a), 



oo 



dr 



Now if we perform the substitution by letting V{r) = V, we have, 

A<p (Vn^) = r' ^2^, where p{V) = -2 ^ [V-\V)]-' (8) 

^ ^ Jv{oo) V^(^) ~^ ^^ 

We can now see that the integrand is linear in p{V). To make the integral look nicer, 
we can let V{x) = Vmax, and applying the condition of asymptotically AdS, we have. 

Equation (9) is in fact a Volterra equation of the First Kind. It is a particular form of 
the generalized Abel equation [17], which is of the form 

yy±-dt = f{x), 0<A<1. (10) 



(x-t) 
and has solution. 



(^\ _ siii(^^) d /•* f{x) dx 



Applying this to equation (9), we have. 



^ 1 rV A(f)[ \/Vmax ) dVmax 

PiV) = -^ / _ (12) 



Now substituting back in for the original variables. 



dr ^ 2 d f^ aA(h(a) , , , 

-^^r =—K7 , ^^ ' da 13 

dV irdVj^ VV-a^ 



1 1 /-V^ aA0(a) 



TT 



, da (14) 

From equation (14), in principle one can solve the integral for any A0(a) and then 
rearrange to find V{r) and thus f{r). 



2.2 Analytic extraction using equation (14): Pure AdS case 

Now that we have an integral inversion technique for extracting /(r), we can test it analyti- 
cally. We will choose the simplest case of pure AdS in global coordinates, i.e. /(r) = l + r^. 

From equation (3) we get A(/){a) = vr in pure AdS, thus using equation (14) we have. 



1 
r 
1 



VW) 



a 



: da 



fir) 



1 A/V"(r) - a^ 

VV{r) - 1 
r^V{r) = 1 + r^ as expected 



(15) 



2.3 Numerical extraction using equation (14) 



Of course not all integrals can be evaluated analytically, but we can solve them numerically, 
which means we can test the extraction of a more complicated f{r). Firstly we can consider 
a spacetime for which the position of the local maximum of V{r) is at r ^ 0: 



0.5 



1.5 2 2.5 3 3.5 



4.5 5 



Figure 2: Extraction of V{r) where /(r) = 1 + r^ 
entire spacetime. 



r2 + l' 



In this case we can extract the 



Now we can try a more complicated /(r) by adding some Gaussian bumps: 
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Original analytic V(r) — 
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Figure 3: Extraction of V{r) where /(r) 



1 + r" 
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A; = (0.05, 0.03, 0.1), fii = (1, 1.3, 2) with a^ = (0.1, 0.1, 0.3). 



j=i 
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In this case we can only reproduce the spacetime down to a mimimum value of r when 
J^' = 0. This is because null trajectories become circular orbits at this point. Null 
geodesies which go past this critical value will end up at the singularity and will thus never 
reach the boundary. 

In summary, we have found an analytic inversion method for extracting static spher- 
ically symmetric bulk spacetimes with one function given information of the bulk-cone 
singularities. Such spacetimes that can be extracted from this method include pure AdS 
and eternal black holes in AdS. The method was found by identifying linearity in the equa- 
tions of motion for null gedesics, and use of Abel's equation. It was shown from graphical 
examples that any spacetime of this type can be extracted down to the local maximum of 
the potential. 

In the next section we will expand on this method to static spherically symmetric space- 
times with two functions. The extraction method will be similar, but we will find that 
information of the bulk-cone singularities will not be enough to determine the metric down 
to the local maximum. 



3 General static spherically symmetric spacetimes 

Metrics of the most general static, spherically symmetric spacetimes in n + 2 dimensions 
are given by. 



ds' = -fir)df + -— + r'dnl (16) 

n(r) 

Applying the condition that the spacetime is asymptotically AdS, we let f{r),h{r) -^ 
1 + r^ as r —>■ oo. 

Following the proceedure in Appendix A, the equations of motion along null trajectories 
are, 



« i _ 1 ^2 _ Kr) 



fif-y r2' /(r 

which leads to the two modified quantities 



'^-^ 



(17) 



At(a) = 2 / 



adr 



Ua) ^f(r) h{r) {a' - ^) 
and 



A0(a) = 2 / ^v"^^^ (19) 



,™H r^^h{r) (a2 _ l^) 
where the same definitions are used as in Section 2. 

3.1 Extracting h{r) 

In parallel with the extraction of section 2, we will consider the A(p{a) function only, only 
now there are two functions to extract to fully determine the metric. We will start by 
considering the extraction of h{r). 

It was noted in the introduction that a relationship has been found between the proper 
length C of boundary-to-boundary spacelike geodesies and the entangement entropy Sa of 
the dual CFT. It will now be shown that given knowledge of S^ and thus C, we can extract 
h{r). 

The expression of the proper length C for boundary-to-boundary static spacelike geodesies 
with the metric given in (16) (see Appendix C) is given by, 

J J ^h{r) (r2 - J2)= 

2r 
If we massage the integral by letting y{r) = ^ we get: 

VHr) 
C{J)= [ -Mk]=dr (21) 



Notice this is a linear integral equation and is in fact a modified version of Abel's 
equation, which has solution, 

, , 2 d r JC(J) , , 
Substituting back in for h{r), we have the final equation with which to extract h{r) given 

' -'' r '^^'^ ij (23) 



A//i(r) r-K dr J ^^ ^ ,p - r^ 

3.1.1 Analytic extraction using equation (23) 

Now we can check the extraction with some trial functions. Since h{r) behaves like 1 + r^ 
for r -^ oo, we begin with the simplest case^^ of h{r) = 1 + r^ . 

We start by determining the proper length C 

£C n = -2 I ^^^ (24) 

Jpo V(l + r2)(r2-j2) 

p'^-'^' dx 



x = r^-J^^C{J) = f 

Jo 



a/eOtTTTj^ 



dJ) = 21n|v5l±^±^IZ? 



Now we can determine h[r), using equation (23) we have. 



(25) 



1 1 d r jcu) , , 1 rf ,, , 

' ^•^ = -x:AO (26) 



^Jh{r) r-K dr J p^^ ^ ,p - r^ ttt (ir 
where 

/(,) = r ^££L in ( VS±T+^ZZ I (27) 

Integrating by parts [18], gives 



\/x — r^ /■'" v' 



X — r 



2 






^Note there is a simpler case with which one can check the inversion, where h(r) = 1, which is Euchdian 
space. In this case C{J) is just the length of a chord on a circle of radius p^, i.e. C{J) = 2-\/pq — J^. Then 
it is just an exercise in integral solving to show that this does indeed imply h(r) = 1. 



where 



' p2 _ r2 



h (r) = -2 \IpI - r2 + 2 Vr2 + 1 arctan | ^ / ^ | (29) 



and 

/2(r) = /2i(r) + /22(r) + /23(r) (30) 

where 



/2i(r) = 2 Jp2 _ r2, /22(r) = -vr Jpg + land/23(r) = TrVr^ + l-2yjr^ + 1 



arcsm 



'pi-r^ 

P^ + 1 
(31) 



J(r) = 2Vr2 + l arctan \ \l ^ T ~ '^'J pI + ^ + ^^^^ + ^ ~ ^^^^ + ^ ^^^^^^ ( \/ ^°2 



TT 



Vr2 + 1 + const. 
1 rf 



— 7rvr2+T 



A//i(r) rixdr 

=^ h{r) = 1 + r^ as expected 

(32) 
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3.1.2 Numerical extraction using equation (23) 

We consider an oscillatory h{r) to show that the extraction method is not limited to simple 
functions, and that the inversion reproduces every bump. 



Extracted h(r) data + 

Original h(r) + 



+ ++4-., ,.-4-^+ 
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Figure 4: Testing equation (23) for h{r) = 1 + r^ — 



cos(2.8r) 



Notice that we can only extract down to h{r) = as this is the position of the event 
horizon. 

3.2 Extracting f{r) 

Now that we have a method of determining h{r), we can use the original method described 
in section 2 using the path of null geodesies. The method will be exactly the same except 
for the modified equations of motion given in (17). 

By substituting in the modified f, using the same substitutions as before, we get from 
Abel's equation [17] that: 






max I ^ ''max 



where q{V) = —2 



dr 



V 



dV^f[V-^V)]h[V-^{V)] 

(33) 
Sustituting the variables back into the equation as before, we finally get: 



V'(r) + 



V(r) 



rIi[V{r)] V h{r) 



0, where h[V] = ---j^ / — J^^ ' (34) 



TT dV Ji /iT" 



a" 
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As we know how to extract h{r), we can solve equation (34) (e.g. by separation of variables) 
to find V{r) and thus f{r). 

3.2.1 Analytic extraction using equation (34) 

Again, using the simplest case of pure AdS, we know that A0(a) = it and h(r) = 1 + r^. 
We can now use (34) to solve for V{r): 



h[V] 



d 



w 



a da 



1 



dV ii y/V-a^ 2^/V^\ 



=> V'{r) 
Separation of variables gives [18], 



2y/V{r){V{r)~] 



In 



2^/V'^{r) - V{r) + 2V{r) - 1 
1 + r^ 



In 



VTT^+i 



VT+r^-i 



V{r) 



/(r) = 1 + r as expected 



(35) 
(36) 



(37) 
(38) 



3.2.2 Numerical extraction using equation (34) 

Here we picked two different functions h{r) and /(r) to fully test the extraction. 



Extracted V(r) data + 
Original analytic V(r) 
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Figure 5: Testing equation (34) for V{r) = 1 -\ — - — 



cos(2.8r) 



1 2 



and h{r) = 1 + r^ 
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As in section 2, we can only extract up to when V{r) is a local maximum. There is no 
limiting in the extracting of /(r) due to extracting h{r) first as the radius of the event 
horizon is smaller than the radius of null circular orbits. 

In summary, we have found an analytic inversion method for extracting static spher- 
ically symmetric bulk spacetimes in three dimensions descibed by two functions, given 
information of the bulk-cone singularities and given the entanglement entropy of the two 
dimensional boundary CFT. Such spacetimes that can be extracted from this method 
include deformations of three dimensional AdS by a self-gravitating gas of radiation (so- 
called "star" geometries). In deriving the method, Abel's equation has been used twice on 
naively independent inversions. This might hint at a relationship between spacelike and 
null boundary-to-boundary geodesies in the context of this extraction. It was shown from 
graphical examples that any spacetime of this type can be extracted down to the local 
maximum of the potential for /(r) and that the extraction of h{r) does not change this 
limiting factor. 
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4 Discussion 

In this paper we have given an analytic method for extracting bulk geometries given the 
location of singular correlators and the entanglement entropy in the boundary CFT in the 
context of the AdS/CFT correspondence. This was achieved in two particular types of 
geometry: Static spherically symmetric, asymptotically AdS spacetimes described by one 
function /(r), and the same type of spacetime described by two independent functions /(r) 
and h{r). This was done by recognising linearity in the expression for A(p{a) (and At{a)), 
therefore allowing an analytic integral inversion. 

One limitatation in the extraction of h{r) was the fact that the area of a minimal surface 
is only the proper length in a three dimensional spacetime. That means our analytic 
extraction given knowledge of the entanglement entropy is limited to the case of three 
dimensions (or a two dimensional boundary CFT). The natural progression would be to 
determine the area of higher dimensional minimal surfaces, and try to extract the metric 
from this, although the extraction method will be a lot more involved than the method 
laid out in this paper. Of course, the method for extracting h{r) is not necessarily limited 
to three dimensions. Since the proper length is a physical quantity, intuitively there should 
be some physical quantity in the boundary field theory which is the dual of this quantity 
in the bulk. If that were found, h{r) could be determined in any number of dimensions via 
the method shown. 

In both cases where the metric is desribed by either one or two functions, the entire 
metric can be recovered if the spacetimes have a local maximum when r ^ 0. On the other 
hand, if the spacetimes admit null circular orbits, then the metric can only be determined 
down to the local maximum of V{r). An extension to this work would be to look at 
cases where the spacetime geometry has less symmetry. This can by done by relaxing the 
spherical symmetry in one or more dimensions and by promoting the functions /(r) and 
h{r) to become time dependent. The ultimate goal would then be to completely extract 
the most general asymptotically AdS spacetime given boundary CFT data. 
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Appendix A: Null Geodesies 

For motion along geodesies, we can just consider the equatorial plane. So when calcu- 
lating the equations of motion, we consider just one angular direction (j). 

The equations of motion for t and (p yeild 2 constants of motion 

i f{r) = E and r^0 = J, \x = — , where A is the afiine parameter 

For null geodesies, ds = 0, thus the equation of motion for r along null trajectories is given 
by 

r^ = E'- ^^i^ (2) 

Dividing by J overall, we can just consider the ratio a = j- to define the null trajectories. 
Thus, we have the 3 equations of motion 

a • 1 n o f(r) , , 

f2 = a2-^ (3) 



/W 



rr> ^ lyA 



To calculate the angular separation, we use the fact that -£: = j and thus 

A0(a) = 2r ^ r^-^ dr (4) 



and similarly 



a 



At(«) = 2 / dr (5) 

where the factor of 2 comes from the fact that the null trajectory goes from the bound- 
ary, down to a mimimun value of r, then back to the boundary at r = oo 



Appendix B: Derivation of equation (4) 

We can calculate the derivative of any function of the form a(a) = J^. s c(a, r) dr, where 
h{a) and c{a, r) are continuous, by using the definition of a derivative, i.e. 

a [a) = lim — ^ -^ 

Sa^o da (-,\ 

poo \^) 

= — h'{a) c{a,h{a)) + I c'{a,r)dr 

Jb(a) 
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This can now be applied to the definitions of A(p{a) and At{a): 

a±At(a) = ^-^^^^^^4^^^^= + 2a T ±(^ ^ \ rJr 

d<^ rlJ«) \/(r„.„(«))V«'-y(r„.„(«)) ir^H da \rW{T) ^ a^ - V[r) ' 

2rV„(«) r- a_ ^^^ 



d 
A(f){a) using eq. (1) 



da 



(2) 



Appendix C: Spacelike Geodesies 

For spacehke geodesies in with metric given by equation (16) we have: 



+2 I ' , ^212 



-fir)t' + —- + r'<j>' = e (1) 

n{rj 

where e > and will be normalised to 1 by redefining the affine parameter A. 

To calcualte the proper length, we can consider the simplest case of static geodesies (i.e. 
t = 0). Hence, 

^2 



h{r] 



+ r^d)'^ = 1 



r =\lh{r) ( 1 ) by subtitution of 






(2) 



The definition of the proper length is given by: 

f f /■'"' dr 

C{J) = ds= dX = 2 — , where f |r-™,„(j)= (3) 

Where r |r™™(j)= =^ rmin{J) = J or h{rmin{J)) = 0. But at h{r) = we have an 
event horizon, and static spacelike geodesies cannot extend inside an event horizon, so we 
use TminiJ) = J- We have also introduced a maximum radius r = po to regulate C. This 
corresponds to the ultra-violet (UV) cut-off in the dual CFT and so is perfectly valid. Thus 
the integration region is po ^ '^ ^ «^- 

Thus the final expression for the proper length is given by 

rPo 2r dr 

^^'^^^'^Jj ^h{r) (r2 - J2y ^^^ 
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